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ABSTRACT 

We present the first definitive measurement of the absolute magnitude of RR Lyrae c-type vari- 
able stars (RRc) determined purely from statistical parallax. We use a sample of 247 RRc selected 
from the All Sky Automated Survey (ASAS) for which high-quality light curves, photometry and 
proper motions are available. We obtain high-resolution echelle spectra for these objects to deter- 
mine radial velocities and abundances as part of the Carnegie RR Lyrae Survey (CARRS). We find 
that My rrc = 0.52 ± O.ff at a mean metallicity of [Fe/H] = —1.59. This is to be compared with 
previous estimates for RRab stars (MyRRab = 0.75 ± 0.13) and the only direct measurement of an 
RRc absolute magnitude (RZ Cephei, My.RRc = 0.27 ± 0.17). We find the bulk velocity of the 
halo to be (WTr,Wg,W z ) = (10.9, 34.9, 7.2) kms -1 in the radial, rotational and vertical directions 
with dispersions (aw* > awo > aw z ) = (154.7, 103.6, 93.8) km s _1 . For the disk, we find (W„, We, W z ) = 
(8.5, 213.2, —22.1) kms" 1 with dispersions (aw„ > aw B > aw z ) = (63.5, 49.6, 51.3) kms~ . Finally, we 
suggest that UCAC2 proper motion errors may be overestimated by about 25%. 
Subject headings: cosmology: distance scale - galaxy: kinematics & dynamics - galaxy: fundamental 
parameters - galaxy: structure - stars: distances - stars: RR Lyrae variables 



1. INTRODUCTION 

Determining distances by use of multiple methods has 
a long and distinguished history in astronomy from antiq- 
uity to the present day. Aristarchus of Samos first deter- 
mined the Moon-Earth distance from the lunar eclipse. 
A century later, Hipparchus checked Aristarchus' val- 
ues using the independent method of terrestrial paral- 
lax: the position of the lunar limb during solar eclipse 
as seen from Alexandria and Hellespont. More recently, 
astronomers have demanded that multiple methods be 
employed for zeroing in on the precise parameters gov- 
erning the currently observ ed and mysterious ac celerated 
expansion of the universe ( | Albrecht et"al]|2006| . 

Pulsating variables have enjoyed a privileged role in 
the local volume. Due to their characteristic light curves 
and relatively bright absolute magnitudes compared to 
the bulk of main sequence stars, they can easily be iden- 
tified and readily measured in nearby galaxies. With a 
local zeropoint for these systems, it is straightforward, 
modulo metallicity and reddening issues, to determine 
the distances to external systems. Cepheids, and their 
more common, albeit fainter, relatives within the insta- 
bility strip, the x (RRL) variables, have been two key 
elements in the historical endeavor to launch humanity 
out of the solar system and Milky Way and into the local 
cosmos. Indeed the Cepheids currently serve as the an- 
chor of the cosmological distance scale, having allowed 
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the most precise measurement the local rate of expan- 
sion of the universe, Hq, to date (Freedman et al. 1994, 
Freedman et al. 2001). In principle, once properly cal- 
ibrated, all stars of known absolute magnitude should 
yield identical measurements of the distances to nearby 
galaxies. This has not historically been the case. 

In particular, the RRL My calibration has varied 
by almost half a magnitude depending on the method 
adopted, and, therefore, consistency with the Cepheid 
distance scale as well as other distance metrics has been 
difficult to establish. As a result, attempts to deter- 
mine RRL absolute magnitudes were largely abandoned 
over the past decade with a few notable exceptions (e.g., 
Dambis et al. 2009). However, recently there have been 
new efforts. Benedict et al. (2011) used HST trigono- 
metric parallaxes of 5 RRL stars (including 4 RRab and 
1 RRc variable) to obtain an average My = 0.45 ± 0.05. 
Klein et al. (2011) recently used mid-IR data from the 
WISE satellite to infer a mid-IR Period-Luminosity re- 
lation. 

In this work, we present a third measurement of RRL 
absolute magnitudes using the method of Statistical Par- 
allax (Sir). The large number of RRL that have been dis- 
covered in the last decade allow us to make a fresh assault 
on this issue. Historically, RRL star distance measure- 
ments from Sir have come in systematically shorter than 
other distance indicators, in particular Cepheids (Barnes 
& Hawley 1986, Hawley et al. 1986, Layden et al. 1996 
(hereafter, L96), Popowski & Gould 1998a, Popowski & 
Gould 1998b , Gould & Popowski 1998 (hereafter, col- 
lectively PG 3 ), Dambis et al. 2009). It is not yet fully 
understood why either this method or these two classes 
of objects should yield different distances to the same 
galaxies. Thanks to automated synoptic all-sky surveys 
like the All Sky Automated Survey (ASAS, Pojmanski 
2002), the number of RRL stars that have reliable light 
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curves has increased by a factor ~5 relative to the pre- 
vious "state-of-the-art" . The ASAS program has identi- 
fied approximately 2000 RRL stars with 300-500 epochs 
of photometry. By obtaining high-resolution spectra for 
these targets, we can both measure the radial velocities 
and metallicities needed for Sir and address outstand- 
ing issues of systematics. The light curves allow us to 
accurately determine pulsation phases and permit the 
measurement of the radial velocity at a single fiducial 
phase at which the pulsation velocity equals the star's 
systemic velocity (see Kollmeier et al. 2009). Tradition- 
ally, obtaining the radial velocity component for RRL 
was laborious, requiring multiple epochs of spectroscopic 
observation. The determination of the phase-velocity re- 
lationship for RRabc variables (Liu 1991, Kollmeier et 
al. 2009, Preston et al. 2011) allows a far more efficient 
strategy for obtaining critical radial velocity information 
with which to compute Sir as we discuss further in Sec- 
tion 3. 

Historically, RRc variables have been either excluded 
from Sit analyses or only approximately analyzed. This 
is primarily due to two factors. First, their hotter 
temperatures make it more challenging to determine 
abundances from low-resolution, low signal-to-noise ra- 
tio (SNR) spectra (see Layden 1994) and, as a result, 
these objects cannot be robustly classified by population 
(halo/disk) as required by modern Sir. However, high- 
resolution echelle observations circumvent this issue and 
allow, for the first time, a definitive S7r analysis from 
RRc variables alone. Second, there are fewer RRcs rela- 
tive to RRabs, and it is only now that samples are large 
enough to perform a robust, self-consistent, pure RRc Sir 
analysis. We analyze our full (RRab + RRc) sample in a 
future work (Kollmeier et al. in preparation) and restrict 
our attention here to our RRc sample. 

In Section 2 we present a brief overview of Sir to remind 
the reader of the basic principles of the technique. In 
Section 3 we present our sample selection, observations, 
data reduction, and analysis methods. In Section 4 we 
review our updated methodology for determining Sir, the 
results of which are discussed in Section 5 and compared 
to previous Sir results in Section 6. Finally, in Section 
7 we discuss our results in light of recent and historical 
works on the absolute magnitude scale of RRL variables. 

2. STATISTICAL PARALLAX 

The basic principle that underlies statistical parallax 
is that the absolute magnitude of any stellar population 
characterized by a particular velocity ellipsoid should 
have a single true value when derived from either trans- 
verse or radial kinematics of that tracer population. The 
radial velocity (RV) determination of the ellipsoid is in- 
dependent of any assumptions about distance, but the 
transverse determination requires an assumed value of 
the absolute magnitude. The RVs alone yield values 
for the the velocity ellipsoid in units of kms -1 . The 
proper motions, when scaled by the square root of the 
flux (^scaled — M x 10 v / 5 ), yield values for the veloc- 
ity ellipsoid in units of masyr -1 . The ratio of these 
two is a distance, D, which is the distance that an 
RRL star would have to be in order to have V = 0. 
And hence, the absolute magnitude of the RRL stars is 
M v = -51og(D/10pc). 



2.1. Basic Formulae 

The basic formulae for computing statistical parallax 
in the presence of observational errors have been well- 
established (Clube & Dawe 1980a,b; Murray 1983; Haw- 
ley et al. 1986; Strugnell, Reid, & Murray 1986, PG 3 ). 
The method involves a 10-parameter maximum likeli- 
hood fit to the kinematic and photometric data: the dis- 
tance scale r], the 3 first moments ("bulk velocity") W{ 
of the sample velocity distribution, and the 6 indepen- 
dent second moments, <r?-. The key parameter of interest 
for distance scale determination is the value of the true 
versus assumed fiducial absolute magnitude of the tracer 
population. 

T) = io( Mfld ™^ _Mt ™V 5 (1) 

The error in determining this distance scaling is given by 
(Popowski & Gould 1998a): 
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and the last expression is for the (more intuitive) 
isotropic case o~ij = crSij, W is the bulk motion, a is 
the velocity dispersion, and 
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where Nrv is the number of stars with accurate RV mea- 
surements and N pm is the number of stars with accurate 
proper motion measurements. The evaluation in Equa- 
tion pi is for k — 2.2, which is typical of halo RRL 
samples. 



2.2. Observational Requirements 
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The distance scale accuracy improves as N eS ' , i.e 
the sample of stars with accurate proper motions and 
radial velocities is increased in size. ASAS has pro- 
vided a public catalog of R RL stars identified via high- 
cadence light curve an alysis ( Szczygi el fc Fabrycky|2007 



Szczygiel et al. 2009). In addition to photometry from 
the ASAS survey, these stars have proper motions from 
the second USNO CCD Astrograph Catalog (UCAC-2; 
Zacharias et al. 2004), which covers the entire southern 
hemisphere. Radial velocity information is not published 
for the majority of the Southern Sky and what is there- 
fore required are spectroscopically determined velocities 
for this large sample of objects that has corresponding 
proper motions and photometry. Before beginning our 
extensive observational program, we evaluated the suit- 
ability of the ASAS sample for measuring Sir considering 
known systematic uncertainties. We discuss each of these 
below. 

2.2.1. Multiple Populations 

The key underlying assumption of Sir is that the RRL 
are a faithful tracer population of a single velocity el- 
lipsoid and do not exhibit poorly mixed kinematics, for 
example, from coherent stellar streams. Clearly, as one 
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goes to very large distance (i.e., the outer halo) this 
assumption breaks down as RRL populations become 
increasingly contaminated by satellite debris on non- 
mixed orbits. Indeed, this consideration alone prevents 
us from extending the survey to extremely faint magni- 
tudes where the distances probed preferentially contain 
such kinematics (e.g. Kollmeier et al. 2009). The distri- 
bution of distances for the ASAS RRc sample is shown 
in Figure [I] As can be seen, the majority of our sample 
lies within a distance of 4kpc. The solar neighborhood 
on these scales is known to be sufficiently smooth that 
our statistics are unaffected by kinematic substructure 
(e.g., Gould 2003). 

2.2.2. UCAC-2 Proper Motions 

The proper motion errors from UCAC-2 are typically 
5-6 masyr -1 for the majority of the sample correspond- 
ing to ~ (26kms _1 ) x (D/kpc). At first sight, the proper 
motion errors may seem too large to be useful, given that 
the median distance is about Z? mc d = 2.7 kpc (see Fig. [I] 
determined using an initial value for the absolute mag- 
nitude scale of My = 0.5), and the direction-averaged 
dispersion is about a ~ 120 kms~ . Together, these im- 
ply that the errors are typically a large fraction of the 
quantity being measured. In fact, S7r is extremely ro- 
bust with respect to observational errors. Equation (17) 
of Popowski & Gould (1998a) (hereafter PG98a) shows 
that the statistical error should increase by only a few 
percent relative to Equation (J2J, given these errors. We 
verify that this is actually the case in Section [5] Equa- 
tion (18) of PG98a shows that one must be careful about 
systematic errors that can result from misestimating the 
size of these errors. Indeed, the Sir methodology contains 
an internal check on the proper motion errors. 




Distance (kpc) 

Figure 1. Distribution of distances for sample stars. 

2.2.3. Reddening & Photometry 

Since Sir uses apparent magnitudes, it is important to 
have robust extinction measurements for objects in the 
sample. We use the reddenings provided by Schlegel, 
Finkbeiner & Davis (SFD, 1998) and a standard extinc- 
tion coefficient, Ry = 3.1 to determine Ay Because 
many of our sources are at high Galactic latitude (see 

6 We note that updates to the SFD extinction maps have been 
provided by Schlafly &; Finkbeiner (2012) and Peek & Graves 



Figure [2b, this directly gives adequate values for these 
sources^We adjust the reported reddening values to ac- 
count for the finite distance from the plane of our ob- 
jects. In some instances, our objects are too close on the 
sky with other field stars in projection to yield reliable 
ASAS photometry (see Pojmanski (2002) for details of 
ASAS photometry). These objects are eliminated from 
our analysis. 

With secure photometry and transverse kinematic 
measurements in hand, we began our spectroscopic sur- 
vey. 

3. OBSERVATIONS & DATA REDUCTION 

The observations presented here were made in 2011 and 
2012 with the echelle spectrograph mounted on the 2.5m 
du Pont telescope at Las Campanas Observatory as part 
of the larger Carnegie RRL Survey (CARRS; Kollmeier 
et al. in preparation). Upon completion, our survey 
will have moderate signal-to-noise ratio spectra for ap- 
proximately 1200 RRL stars observable from the South- 
ern Hemisphere. The sky distribution of our full survey 
sample is shown in Figure [21 The grey points show RRab 
stars and the blue points show RRc stars. Observations 
were designed to reach SNR^ 15 in the order containing 
the Mg I triplet at 5170A at the target phase. Exposures 
for a given target bracketed observations of a ThAr lamp 
through a 1.5x4-arcsec slit. 




Figure 2. Sky distribution of our ASAS targets. The RRab sam- 
ple is shown in grey points and the RRc sample sample is shown 
in blue. The analysis presented in this paper includes 247 of the 
RRc sample points. 



3.1. Synoptic Observations 

The photospheric velocity of RRL stars changes by 
many kms -1 over the course of a pulsation cycle. Tra- 
ditionally, one has taken multiple observations through 
the pulsation cycle and fit the resultant velocity func- 
tion (when the phase is known) in order to obtain the 
systemic velocity of the target (e.g. Layden 1994). We 
adopt the time-averaged velocity of the pulsation curve 

(2011). These new maps qualitatively agree, sharing the conclu- 
sion that the SFD extinction map is generally thought to be over- 
estimated, but do not quantitatively agree in terms of the precise 
magnitude of this difference (14% and 2% in each study). The ef- 
fect of adopting the more extreme correction leads to a correction 
in 77 of A?7 = 0.03 x log(10)/5 = 0.014. We adopt the SFD values 
until a definitive calibration is agreed upon. 
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as the cen ter-of -mass ve l ocity of the star as detailed in 
Liu (1991) and Preston (2011) for RRabs. Integration 
of detailed velocity curves ol the RRc variables T Sex, 
TV Boo, DH Peg, and YZ Cap, all measured with the du 
Pont echelle, shows that the pulsation velocity is equal 
to the star's time-average velocity at phase 0.32, reck- 
oned relative to maximum light. Our observations were 
all made as close to this phase as possible, and velocity 
corrections were applied adopting the following correc- 
tion: 

Acf> RRc = 62.0kms- 1 (0.32 - <f> obs ), (5) 

Figure [3] shows the resultant phase distribution for our 
targets. Note that these velocity corrections are no larger 
than 5kms _1 for any given star. 




Figure 3. Phase distribution of our sample of 247 RRc variables. 
The red line in the figure shows the phase targeted for the RRc 
sample throughout our survey. 



3.2. Data Reduction Pipeline 

In order to ensure uniform quality across our sizeable 
database, a data reduction pipeline was constructed fol- 
lowing Kelson (2003). The pipeline divides the reduction 
process into first processing the calibration frames, and 
then using these to reduce the science frames and extract 
spectra. The calibration stages involve processing the 
frames for flat-fielding, determining the order edges and 
y-distortion using sky flats, and obtaining wavelength so- 
lutions using a Th-Ar lamp frame. Once these frames are 
produced, the science spectra are reduced. The science 
spectra are first divided by the processed and normalized 
flat-field and then wavelength calibrated. The spectra 
are then sky-subtracted and the spatial profiles of the 
resultant 2-D spectra are used to find extraction aper- 
tures. The pipeline takes advantage of parallel process- 
ing, which speeds up the reduction process significantly 
for a large number of objects. 



Post-extraction processing of the spectra was done 
with the IRAlQ ECHELLE package. Radial velocities 
were measured relative to a high SNR template of the 
metal-poor star HD140283 using the IRAF FXCOR rou- 
tine. The template spectrum was taken with the du Pont 
Echelle spectrograph in the same configuration as survey 
stars. The heliocentric radial velocity for HD140283 has 
been determined, independently in four separate high- 
resolution studies, to be Vheiio = —170.92 ± 0.19 kms -1 
(e.g. Tsangarides et al. 2004, Aoki et al. 2002, Lucatello 
et al. 2005, Latham et al. 2002). The cross correlations 
were made on the wavelength interval 4900A - 5500A. 
Objects that were originally classified as RRL stars based 
on their light curves but, upon inspection of their spec- 
tra, were found to be binaries or other non-RRL variables 
were removed from the sample. 

The radial velocity error estimates returned by fxcor 
are extremely small, so we estimate the true errors by 
making repeat measurements of the corrected velocities 
on a subset of 11 stars. From these we determine ctrv = 
1.96 km s -1 . We add this in quadrature to our mostly 
tiny formal errors. We note that, from PG3, the velocity 
errors enter the final result as Arj/rj ~ (c T Rv/o'haio) 2 /6 < 
10~ 4 and so have no practical effect in any case. We 
nevertheless include this correction for completeness. 



3.3. Abundance Determination 

The measurement of absolute abundances for large 
samples of RRc stars has not been undertaken previously. 
We therefore must calibrate our pipeline abundance mea- 
surements (which yield consistent relative abundance 
determinations) to the scant data currently available. 
Based on analysis of the high SNR du Pont echelle spec- 
trum of RRc variable YZ Cap, we adopted a single set 
of atmospheric parameters (T c ff =7000K, \ogg = 2.2, 
^micro = 2.5kms~ 1 , [a/Fe]=+0.35) for all survey stars. 
The value of T e g is somewhat constrained by the small 
area occupied by RRc stars in the color-magnitude dia- 
gram, and our survey spectra are not of sufficiently high 
SNR to determine this parameter independently. The 
abundance measurements in RRL stars are relatively in- 
sensitive to the adopted values for logg and v m - lcro . For 
our final calibration, we rely primarily on a detailed 
study (Govea et al. in preparation) that examines the 
sensitivity of RRc abundance determinations with re- 
spect to stellar parameters using a small sample of high 
SNR spectra also obtained with the du Pont echelle. This 
is the only study we know of in which the sensitivity of 
RRc abundance measurements are systematically stud- 
ied. We also compare to the known distribution of RRab 
metallicities, as a secondary calibration. We defer a more 
detailed analysis of this metallicity calibration to a fu- 
ture work where more high SNR, high-resolution data 
are available for comparison (Kollmeier et al. in prepa- 
ration). 

Pipeline abundance measurements were performed by 
fitting to a grid of synthetic spectra generated by 

7 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under a cooperative agreement 
with the NSF. 
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MOOGH Spectra reduced by our spectroscopic pipeline 
were cleaned of remaining noise spikes and smoothed 
using a 3-pixel boxcar filter. Owing to the relatively 
broad lines in RRL stars, this smoothing has essentially 
no effect on the abundance determination. RRL stars 
are warmer than the Sun and generally metal-poor and 
therefore contain few spectral features beyond 5500A. 
As our SNR is maximized between approximately 4400A 
and 5500A, we therefore used two broad spectral ranges 
for our synthetic fits: 4400-4680A and 5150-5450A. The 
latter region contains primarily neutral-species transi- 
tions as well as the Mg I b lines. The former region has 
numerous singly-ionized species including Ba II 4554A. 
The synthetic spectral grid covers a range of effective 
temperature (Teg), surface gravity (logp), metal abun- 
dance and microturbulent velocity (v m i cro ). The atomic 
line lists for the syntheses were begun with the Ku- 
rucz line database, then refined until a good match was 
achieved for the spectrum of the Sun. The smoothed syn- 
thetic spectra were compared with the observed spectra 
and the parameters of those with lowest x 2 were chosen 
as fits. In Figure [4] we show a histogram of our pipeline- 
derived and calibrated abundances averaging the (inde- 
pendently) derived abundances for the two broad spec- 
tral windows. 
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Figure 4. [Fc/H] determinations for our sample of 247 RRc vari- 
ables. Abundances were determined by averaging two broad spec- 
tral windows which had maximal signal to noise as well as numerous 
spectral features. 

4. DISK/HALO SEPARATION 

With the kinematic and chemical information mea- 
sured as described above, we are prepared to segregate 

8 MOOG is a publicly available code to determine abun- 
dances in stars through LTE analysis. MOOG is available at 
http:/ /www. utexas.edu/~chris/moog. html 



our sample into the Halo and Disk populations. We con- 
vert our proper motions and radial velocities to radial, ro- 
tational, and vertical velocities (U, V, W) as per L96. We 
correct these velocities for the dynamical solar motion 
(+9, +250, +7 kms -1 ) and we rotate the Sun-centric 
(U, V, W) 3-space velocity to the local frame of the star 
assuming cylindrical symmetry of the Galaxy. The cu- 
mulative distribution of these rotation angles is shown 
in Figure [5j Note that they span the range of roughly 
±25°. The resultant velocities {v n ,vg,v z ) for our sam- 
ple are shown in the left panels of Figure [6] for direct 
comparison with Figure 3 of L96 which shows the same 
quantities for their sample of RRab stars. Similarly, the 
right hand panel of Figure [6] shows the rotational veloc- 
ity component as a function of measured abundance for 
direct comparison with Figure 4 of L96. It is important 
to remind the reader again that we are using RRc stars 
rather than RRab stars as tracers of the velocity ellip- 
soid, in contrast to L96. Despite this difference in tracer 
population, the kinematics of this sample are similar to 
those of L96. 
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Figure 5. Distribution of rotation angles between the Sun-centric 
and star-centric frames of reference. 

While it is not possible to know with certainty whether 
any particular star is a member of the disk or halo, we 
can clearly see a track of stars that has chemical and 
kinematic properties similar to a disk-like population. 
These stars are relatively high metallicity and are ro- 
tating in the direction of the Sun with similar speed. 
The exact demarcation of the disk/halo line is necessar- 
ily uncertain. Following L96, we therefore perform our 
analysis on subsamples with the aim of determining how 
sensitive our results are to these uncertainties. As we 
will show, they are not. We identify two halo/disk sub- 
samples. The first (HALO-l/DISK-1) considers all stars 
rightward (upward) of the line vg — 400 [Fc/H] — 225 to 
be members of the disk/thick disk and they are excluded 
from the halo population of interest. Of our 247 objects, 
this procedure assigns 34 to the disk and 213 to the halo. 
These objects are shown by the red squares in Figure [6j 
The second (HALO-2/DISK-2) considers all stars with 
[Fe/H]> -1.0 and vg > 100 km s -1 to be disk stars. Ob- 
jects satisfying these criteria are shown as blue triangles 
in Figure M Of our 247 objects, this classification elim- 
inates 30 leaving 217 for our Sn analysis. We note that 
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our HALO-l/DISK-1 and HALO-2/DISK-2 designation 
is slightly modified from L96 given our sample. However, 
our results are relatively insensitive to these distinctions, 
although our errors will necessarily scale with the square 
root of the number of retained objects. 
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Figure 6. Kinematics and abundance information for our 247 
RRc stars. The left panels show the values of v n ,vg, and v z , that 
is, the velocities in the radial direction, the direction of Galactic ro- 
tation, and toward the Galactic pole respectively. These velocities 
are corrected for rotation w.r.t. U, V, W velocities. The right panel 
shows the rotational velocity component as a function of [Fe/H]. In 
each panel, red squares denote the "DISK-1" population and blue 
triangles show the "DISK-2" population. These stars are elimi- 
nated from our ML analysis as described in the text. 



5. RESULTS 
5.1. Halo Velocity Ellipsoid 

We carry out a 10-parameter Markov Chain Monte 
Carlo (MCMC) maximum-likelihood fit to the data. 
These include the nine kinematic parameters describ- 
ing the velocity ellipsoid (three diagonal elements cor- 
responding to the bulk motions and 6 dispersions) and 
one parameter describing the absolute magnitude scal- 
ing rj. Table [I] gives the values and errors of the 10 
Sir parameters ior our samples. For reasons explained 
in Appendix A, their errors are almost perfectly Gaus- 
sian. Therefore the 10-dimensional likelihood surface is 
completely specified by these numbers, together with the 
correlation coefficients, which are also given in Appendix 
A. The fact that the parameters derived from these two 
different fits are consistent at well below 1 a shows that 
our choice of disk/halo separation does not significantly 
affect our results. 

For our HALO-1 sample, we obtain a value of the ab- 
solute magnitude of RRcs of My = 0.47 ± 0.10 and for 
our HALO-2 sample, this value is M v = 0.49 ±0.10. We 
discuss these values in the context of previous measure- 
ments in Section [6j 

5.2. Test of Accuracy of UCAC-2 Proper-Motion Error 

Estimates 

In the analysis reported above, we adopted the proper- 
motion error estimates in the UCAC-2 catalog. As dis- 
cussed by PG98a, if the size of these error estimates were 



systematically too big, it would induce a systematic er- 
ror in our distance-scale estimate, making the RRL ap- 
pear to be systematically further (hence more luminous) 
than they actually are. This is a potential concern be- 
cause Zacharias et al. (2004) report external tests only 
on their mean proper motion estimates and not on their 
error estimates. Fortunately, we are able to perform an 
internal test on these error estimates. This test shows 
first that the UCAC-2 error estimates are probably over- 
estimated by about 25%, and second that this level of 
mis-estimation has almost no effect on our results. We 
note that our RV errors are too small to have any effect, 
whether correctly estimated or not. They will therefore 
be ignored in the following discussion. 

We first review the basic physics that permit such a 
test and then present results. What is now called "sta- 
tistical parallax" was formerly (in the first half of the 
last century) divided into two effects. One, called "sec- 
ular parallax", compared the mean bulk motion (Wi) in 
RV with that in (flux-scaled) proper motions. The other, 
also called "statistical parallax" compared the amplitude 
of the dispersions (o~ij) in RV and proper motions. In 
modern S7r, these are done simultaneously in a single fit. 
However, the errors in the proper-motion error estimates 
enter very differently into these two components. For 
the aij comparison an overestimate of the errors leads 
to a corresponding underestimate of the intrinsic proper 
motions of the stars, which can be reconciled with the 
dispersions measured in RV only by placing them further 
away. However, for the Wi measurement, which is first- 
order (not second-order) in the proper motions, there is 
no such effect. If one wished to be completely safe from 
any such error-bar misestimates, one could in principle 
choose to just compute the "secular parallax" . 

However, the maximum likelihood formulation of Sir 
permits a more sophisticated approach. The likelihood 
function is a direct and sensitive indicator of the consis- 
tency of the two components of the Sir determination 
when the error bars are systematically rescaled. Ta- 
ble [2] shows results for several such rescalings. The log- 
likelihood is maximized for a rescaling factor / err = 0.5, 
while f CII = 1 yields a log-likelihood that is lower than 
this AlogL = 3.8. The probability of this occurring 
by chance is about exp(— AlogL) ~ 2%. While this is 
certainly not an impossible statistical fluctuation, nei- 
ther is there any compelling prior that the UCAC-2 er- 
rors are exactly as given. We note that at / orr = 0.8, 
AlogL = 1.5, which is very plausibly due to statistical 
fluctuations. Because we do have some prior that these 
errors are not grossly overestimated, and because there 
is no statistical evidence that / err < 0.8, we adopt this 
as our best estimate of the errors. We note from Table 1 
that the result of changing / err from 1.0 to 0.8 is to reduce 
the distance scale by approximately 1%, which is about 
20% of our statistical error. The remaining parameter 
estimates barely change as well. 

5.3. Small Correction for Malmquist Bias 

Because the sample is magnitude limited and the target 
population has an intrinsic dispersion in absolute magni- 
tude, o~m v 1 the sample contains more stars with brighter- 
than-average luminosity than fainter-than-average. This 
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Table 1 

Resultant Parameters and Errors. 



Sample 


V 
(er n ) 


W.l 

(o- TO1 ) 


U>2 
(c™ 2 ) 


IU3 

(°w 3 ) 


(cr 1/2) 


r 1 / 2 
t 22 ^ 

°22 


r 1 / 2 

°33 


Cl2 


Cl 3 


C23 



HALO-1 1.007 11.8 32.5 8.3 155.6 102.4 93.8 0.037 -0.006 -0.114 

0.048 11.1 9.9 6.9 10.2 6.0 5.6 0.075 0.075 0.076 

Anaiytic-c„ 0.045 10.7 9.6 6.4 8.9 5.8 5.3 0.069 0.069 0.069 

HALO-2 1.002 10.1 37.4 6.1 153.8 104.8 93.8 0.018 0.008 -0.144 

0.048 10.9 10.0 6.8 10.1 6.1 5.6 0.074 0.075 0.075 

DISK-1 0.728 4.0 217.1 -24.9 62.1 42.6 49.4 0.000 -0.192 0.008 

0.142 11.3 9.5 10.8 10.2 8.7 12.8 0.198 0.193 0.207 

DISK-2 0.796 13.1 209.3 -19.2 64.9 56.5 53.1 -0.020 -0.397 -0.088 

0.158 12.6 12.9 11.7 11.2 12.2 14.1 0.203 0.190 0.207 



Table 2 

Effect of Proper Motion Error Rcscaling 



Scale Factor Alog£ r\ w\ u>2 w% cj| 2 ^33 2 C12 C13 C23 



1.0 


3.91 


1.007 


11.8 


32.5 


8.3 


155.6 


102.4 


93.8 


0.037 


-0.006 


-0.114 


0.8 


1.51 


0.994 


12.0 


34.1 


8.2 


155.8 


102.9 


94.4 


0.033 


-0.003 


-0.106 


0.7 


0.74 


0.987 


12.1 


35.0 


8.2 


155.8 


103.1 


94.8 


0.031 


-0.002 


-0.102 


0.6 


0.30 


0.981 


12.2 


35.7 


8.2 


156.0 


103.3 


95.4 


0.029 


-0.001 


-0.098 


0.5 


0.09 


0.974 


12.4 


36.6 


8.1 


155.9 


103.4 


95.8 


0.028 


0.000 


-0.094 



leads to a correction 



Arj/rj 
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(6) 



We modify the PG98a estimate to obtain <tm v =0.1 as 
follows. PG98a began with the observed dispersion in the 
Large Magellanic Cloud (LMC) of 0.17 mag. They esti- 
mated that 0.09 mag was due geometric dispersion, and 
so estimated an intrinsic dispersion of 0.14 mag. We fur- 
ther note that for an estimated metallicity dispersion of 
0.5 mag, and assuming a slope of K = 0.214 mag dex -1 
in the RRL My-[Fe/H] relation, an additional 0.1 mag 
of scatter can be accounted for by metallicity variation. 
Subtracting this in quadrature, we obtain 0.1 mag intrin- 
sic dispersion. This leads to a correction A77/77 = —0.006. 

Note that PG98a identified another effect due to the 
dispersion of absolute magnitudes, which scales Arj/rj ~ 
Q.2hcr\ Iv ~ 0.003. We briefly describe this effect and why 
it is proper not to include it here. If stars do not have 
exactly the same absolute magnitude, then their mean 
square scaled (by flux) proper motions will be greater 
than their squared scaled (by mean flux) proper motions. 
This will mimic a larger proper-motion dispersion and so 
cause the stars to appear closer (so dimmer) than they 
actually are. However, because the mechanism of this 
effect is similar to that cause d by misestimation of the 
proper-motion errors (Section 5.2), this small correction 
is already absorbed into the correction adopted there. 

5.4. Final Result 

To obtain our final result, we incorporate the small 
corrections for Malmquist bias and for probable overes- 
timation of the UCAC-2 proper motion errors and take 



the average between our two halo samples, 

M v = 0.522±0.106±0.031 at ([Fe/H]) = -1.59 (7) 

The systematic error (0.031 mag) is determined by com- 
bining in quadrature the variance from different defini- 
tions of the halo (0.010 mag), uncertainty in the degree of 
overestimation of the proper motion errors (0.025 mag), 
and uncertainty in the extinction scale (0.014 mag). 
There is also a systematic error of about 0.05 dex in the 
metallicity scale, and so the mean metallicity at which 
this estimate is valid. 

Because the statistical errors are about 3.5 times larger 
than the systematic errors, a detailed investigation of the 
latter is not warranted at the present time. This issue 
will become somewhat more pressing when we combine 
the analysis of CARRS RRab and RRc stars. 

Table 1 shows a comparison between the errors derived 
from the the MCMC and those predicted analytically 
using the formulae of PG98a, which assume zero mea- 
surement errors and isotropic sky coverage. The actual 
errors are only slightly larger than the ideal errors despite 
the relatively large proper motion errors. In Appendix 
A, we derive the corresponding analytic estimates of the 
correlation coefficients among the 10 Sit parameters and 
show that these are in good agreement with the MCMC 
determinations. 

5.5. Disk Kinematics 

We can run our maximum likelihood machinery just 
as easily on a pure disk population as we can on a 
pure halo population, albeit with the substantially re- 
duced numbers in our disk sample. We perform this 
analysis and include the results in Table [T] Interest- 
ingly, the value we obtain for absolute magnitude (av- 
eraged between the two "disk" definitions) is My = 
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1.08 ± 0.43 at [Fe/H]= -0.7. After taking account of 
the slope K = 0.214 mag dex -1 in the My- [Fe/H] rela- 
tion, this is consistent at the 1 a level with our result 
for the halo sample. Also of interest are the detailed 
values of the velocity ellipsoid for the disk population: 
(W v ,W e ,W z ) = (8.5, 213.2, -22.1) kms -1 with disper- 
sions (crw w > cr We > u Wz ) — (63.5, 49.6, 51.3) kms - . In- 
deed, it is thick-disklike, which gives us confidence that 
we are in fact removing a kinematically distinct popula- 
tion via our chemodynamical criteria. 

6. COMPARISON WITH PREVIOUS PARALLAX 
ESTIMATES 

6.1. Absolute Magnitude of RRL 

Since the work of Clube & Dawe (1980a, b) there 
have been multiple attempts to derive the RRL abso- 
lute magnitude scale from Sn (e.g. Hawley et al. (1986), 
Strugnell, Reid & Murray (1986), Layden et al. (1996), 
PG 3 , Luri et al. (1998), Dambis (2009)). The most re- 
cent Sir estimate of My^^ (and so directly comparable 
to the work presented here) is by PG 3 who found 

My =0.75 ±0.13 at ([Fe/H]) = -1.6 (8) 

from a sample of 182 halo RRab stars taken from Layden 
et al. (1996J3 

Benedict et al. (2011) used HST to obtain trigono- 
metric parallaxes for 5 RRL (4 RRab and 1 RRc). 
Combining their reported values and measurement er- 
rors (together with their adopted intrinsic dispersion of 
o~m v = 0.0577) yields 

My = 0.443 ± 0.067 at ([Fe/H]) = -1.50 (9) 

while, using o~m v = 0.1 (adopted here) yields My = 
0.426 ± 0.080 at ([Fe/H]) = -1.52. To' obtain these 
results, we weight the individual measured magnitude 
offsets from the mean by (o- 2 Mv + ^a?^,^) -1 , where 

K = 0.214 is the slope of the My- [Fe/H] relation. Note 
that Benedict et al. (2011) incorrectly quote somewhat 
smaller errors. See Appendix B. 

Efforts to measure My from RRc stars have been much 
more limited. Hawley et al. (1986) and Strugnell, Rcid 
& Murray (1986) obtained Sir results from 17 and 26 can- 
didate RRc stars respectively, but were forced to make 
a restricted analysis because of small number statistics, 
which they reported only for "completeness". In addi- 
tion, as mentioned above, Benedict et al. (2011) ob- 



tained My = 0.27 ± 0.17 for a single RRc, RZ Cep at 
[Fe/H]=-1.77. 

6.2. Velocity Ellipsoid 

Our results on the kinematics of the halo 
{{W v ,We,W z ) = (10.9, 34.9, 7.2) kms" 1 with dis- 
persions (ctvv„ , o~w e j &w z ) — (154.7, 103.6, 93.8) kms -1 ) 
are in excellent agreement with previous work. Compar- 
ing with PG 3 , our results are consistent at well under 
1-cr despite an entirely different sample. Comparing 
with the most recent (and largest) Sir analysis of 
Dambis et al. (2009) our results are also generally in 
excellent agreement, as almost all parameters of the 
halo and disk velocity ellipsoids agree within 1 or 2-cr 
with the exception of vertical velocity dispersions of 
the halo which are in tension at the 3-cr level. We do 
not know the origin of this discrepancy, but will have 
more leverage to investigate this with our larger RRab 
sample. 

7. DISCUSSION AND CONCLUSIONS 

We have performed the first decisive analysis of the 
absolute magnitude for RRc variables via statistical par- 
allax using the first data from the Carnegie RR Lyrae 
Survey (CARRS). Our current measurements for RRc 
variables yield a 5% distance error which is similar to 
that obtained from modern techniques applied to RRab 
samples. We find a velocity ellipsoid for our disk and 
halo population that is in good agreement with previous 
measurements . 

Already, CARRS provides competitive distance accu- 
racy to other surveys and techniques. At the conclu- 
sion of CARRS, we anticipate a factor of ~ 4 increase 
in the number of tracers and, consequently, 2% distance 
errors. In future work, we will analyze this far larger 
database and have the statistical potency to divide our 
sample into finer mctallicity and kinematic bins than can 
be done presently. This will allow a precision measure- 
ment for comparison with other techniques with the hope 
of a "unified" RRL distance scale. 

In the Gaia era, where space-based parallaxes will 
be available for many of these objects, our database 
of high-resolution spectra should provide useful com- 
plementary information for going beyond distances and 
gaining further understanding of RRL as astrophysical 
objects, rather than merely as test particles. 



APPENDIX 



A: ANALYTIC ESTIMATES OF Stt CORRELATION COEFFICIENTS 

In the main text, we showed that the analytic error estimates derived by PG98a for the 10 Sir parameters closely 
approximate the numerical errors in our sample, despite the fact that the analytic treatment assumes zero measurement 
errors, uniform sky coverage, and no rotation with Galactic position. Here, we extend the PG98a analytic treatment 
to the off-diagonal elements of the parameter determinations and compare the resulting correlation coefficients to our 
numerical values. 

Integrating Equations (30)- (35) from PG98a over N objects uniformly distributed on the sky, we obtain the inverse 



9 Dambis et al. (2009) computed the Sn for a sample of 364 
Galactic RRL stars in the Northern Hemisphere from targets in 
the Beers (2000) catalog using 2MASS photometry to correct the 
infrared-inferred period-luminosity relation for RRL stars, finding 



(M Kg (Adopted)} = -2.33 log(Pp) - 0.818± 0.081. This work does 
not estimate My rr and can thus not be directly compared here. 

10 Luri et al. (1998) found a similar value from a smaller sample 
(144 stars total) than PC 3 . 
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covariance matrix of the 10 Sn parameters ao . . .ag = (In t], Wi, Cu, Ci^j) 

Ko = b 0u = NA V ; b mn = NB m 8 mn (Al) 

where 

A =4+— . A 1 _3 = -tt ? 7 1 , A 4 - 6 = - — , A 7 _ 9 = 0, (A2) 

12 1 

B i-3 = 7^, -^4-6 = 7^2 , -67-9 = ^ ^ 7 (A3) 

and where we have adopted a reference frame in which CV,- is diagonal. (See Equations (37)-(42) of PG98a.) 
Matrices of this form can be analytically inverted, c = b^ 1 as 



coo 



. . 42 \ -1 A X A A 

a o^Z^~r~! ' c 0m = c m0 = - — Cqo, c mn = — h c 00 . (A4) 



PG98a (Equations (43)-(46)) have already evaluated the diagonal elements of this matrix (i.e., the variances of the 
parameters) 

„ , a var(Wi) 1 / 4 Wf \ var(C ll ) 1 / 8 \ var(CU 7 ) 1 , 4 . 



where a 1 = 4/3 + (2/9)k 2 . Here we use Equation) A4) to evaluate the off-diagonal elements, or equivalently the 
correlation coefficients 

cc(lnf7, Wi) = Q w ,i cc(ln?7, Cu) = Q c ,i (A6) 



cc(Wi, Wj) = QwaQw, 3 cc(Wi, Cjj) = Q w ,iQcj cc(C ii; C i3 ) = Q c ,iQc,j (A7) 

where 

/ 9(7- \ ~ 1/2 / 9 \ ~ 1/2 

g^^sgn^^l + ^j , ^-(l + ^J . (A8) 

and all other terms vanish. Note that in making these evaluations, one must return to the "Sun frame" , i.e. add 
(9, 250, 7) kms^ 1 back to the W*. 

In the matrix below, we compare the analytic correlation coefficients (above diagonal) to the actual ones (below 
diagonal). There is good overall agreement. 





1.00 


0.01 


-0.68 


0.01 


0.53 


0.53 


0.53 


0.00 


0.00 


0.00 \ 
0.00 


-0.02 


1.00 


-0.01 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 




-0.66 


0.04 


1.00 


0.00 


-0.36 


-0.36 


-0.36 


0.00 


0.00 


0.00 




-0.04 


0.00 


-0.04 


1.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 




0.60 


-0.01 


-0.40 


-0.03 


1.00 


0.28 


0.28 


0.00 


0.00 


0.00 




0.38 


-0.01 


-0.25 


0.00 


0.23 


1.00 


0.28 


0.00 


0.00 


0.00 




0.41 


0.00 


-0.27 


-0.02 


0.25 


0.15 


1.00 


0.00 


0.00 


0.00 




-0.16 


0.01 


0.11 


0.01 


-0.08 


-0.06 


-0.08 


1.00 


0.00 


0.00 


V 


0.00 


-0.01 


-0.01 


0.00 


0.00 


0.01 


-0.01 


-0.10 


1.00 


0.00 


0.08 


0.00 


-0.06 


-0.02 


0.06 


0.00 


-0.02 


0.01 


0.03 1.00 / 



We note that the errors in these parameters are almost perfectly Gaussian, and therefore completely described by 
their means and covariances. This can be seen as follows. In any linear fit, the parameter estimates, Oj can be written 
as a linear function of the data yv. <ij = Cijdj, dj = J2ki &klfj(tk)yi> where fj{tu) are the trial functions, Bm is the 
inverse covariance matrix of the data, c = b~ 1 , and fry = J2ki Bklfi(tk)fj(tl)- Therefore if the data yi have Gaussian 
errors, then so do the parameters et^. Note that this statement does not depend in any way on the central limit 
theorem, as is sometimes supposed, but only on the fact that linear combinations of Gaussians are Gaussians. 

In the present problem, — 2hiL looks similar in structure to a standard linear-fit x 2 except, first, the parameters 
appear in non-linear combinations, second, the second moments of the velocity ellipsoid appear in an additional term 
containing a log determinant, and third, the scaling factor 77 also appears in a log term. Nevertheless, the linearized 
problem (treated explicitly by PG98a) is extremely close in structure to a standard linear-fit x 2 ■ Therefore, one expects 
similar mathematical properties. This is the underlying reason that the PG98a linearized analysis matches numerical 
results so closely. 
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B. APPROPRIATE WEIGHTING IN DERIVING UNCERTAINTIES 

Here we derive the appropriate weighting scheme for an ensemble of RRL parallax measurements with different 
errors in both My and [Fe/H]. 

Let us consider n stars with measured absolute magnitudes Myi and errors Oi and each with perfectly known 
metallicity [Fe/H]^. And let us initially assume that the slope of the [Fe/H] — My relation is known but the zero-point 
is not. That is, 

My, pred ([Fe/H]) = a + if ([Fe/H] - Q) (Bl) 
where Q is some arbitrarily chosen fiducial metallicity, if is the known slope and a is the unknown zero point. Then 

A /M v , pred ([Fe/H]J - M v ^ x " 



* 2 = E( 

i=l v 

This x 2 is minimized by setting its derivative to zero, i.e., 



aJ^Wi =J2MMv,i ~ ^([Fe/H] 4 - Q)), 



Wi = a, 



(B2) 



(B3) 



We now choose the "arbitrary" fiducial metallicity so that the relation is independent of the choice of if. This can 
be done if we choose _ 

Q = S^M, (B4) 



in which case \ 2 is minimized at 



a = 



E* w l M Vi 

E l w i 



regardless of if. The error in this estimate is the point at which % 2 rises by unity, i.e. 

-1/2 



a (a 



d 2 x 2 /da 2 



(E 



(B5) 



(B6) 



Hence, Equation(B4) gives the effective metallicity at which the measurement is made. 

Now let us suppose that the metallicities are given to us with error bars Sj, so that we can consider simultaneously 
fitting for both the zero-point (now called do) and the five true metallicities, called dj. We can write x 2 



x 2 = 



E 



Oi - [Fc/H], 



Setting the n + 1 derivatives of \ 2 to zero (wrt to do, cij) yields the equation 

TO 



where 



d = MMv,i + KQ), di = Wi{M v ,i + KQ)K + W^Fe/H^ W t = 5T 2 
boo = A b 0i = bio = Ai bij = SijBi 



and 



a^E 



w, 



Ai = Kw z 



Bi = K 2 u h + Wi 



(B7) 

(B8) 

(B9) 
(BIO) 

(BH) 



8=1 



The inverse of this matrix, c = b 1 is given by Equation) A4 1 , which then allows us to evaluate a Q = E„ c 0vd v , 

bo 



a o — coo 



Wi(M v ,i + KQ) ( 1 - ^if ) - Wi [Fe/H] 



This simplfies to 



A M v<i - if ([Fe/H] i - Q) 
fl0 ^ C00 ^ (if E,) 2 + of 



i=l 



^o 1 — E 



(B12) 



(B13) 
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where we note that c no is the error in ao (e.g., Gould 2003a). This looks identical to our previous expression that 
ignored the [Fe/H] errors, except the inverse weights are increased fractionally by (KY>/<r) 2 . This result formally 
confirms one's naive idea that the metallicity error is "equivalent" to an additional error in the absolute magnitude, 
propagated by the slope (K) of the relation. 

Note, however, that in contrast to the case of perfectly known metallicities, one cannot enforce complete independence 
of the result from choice of K simply by adopting Q as the average metallicity weighted by [(KE) 2 + er 2 ] -1 , since the 
weighting itself depends on K. However, in the present case (ICE /a) 2 <C 1, so this has no practical impact. 
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